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Expressions (in closed form) are derived for the first five derivatives of the Gruneisen and Einstein 
functions. Recursion formulas for the successive derivatives of both functions are also given. Computa- 
tions are made and tabular results are given for the first five derivatives of the Gruneisen and Einstein 
functions. 
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1. Introduction 

According to the solid state theory of matter, the mathematical expression for the electrical 
resistivity p, of good conductors is (ref. [1 1 1 ) 



AT 



(20 



(1) 



where, 



A = constant 
# = Debye temperature 
T= absolute temperature 

The function G°(0lT) is frequently referred to as the Gruneisen function. This function is 
approximately unity for T > 6, which indicates that at high temperatures electrical resistivity in- 
creases linearly with temperature. At very low temperatures where TO, the Gruneisen function 
can be expressed approximately by G (OIT) ~ B(T/0) 4 , which indicates that at very low tempera- 
tures electrical resistivity increases with the fifth power of temperature. 

For various calculations, the successive derivatives of the electrical resistivity function may 
be needed. This requires a knowledge of the derivatives of the Gruneisen function. Also, it may 
be possible to make further refinements in the electrical resistivity expression by expanding the 
resistivity function in the derivatives of the Gruneisen function. It was observed that the derivatives 
of the Gruneisen function contain the Einstein function and its successive derivatives. The objective 
of this writing is to obtain expressions for the successive derivatives of the Gruneisen and the 
Einstein functions. 



'Figures in brackets indicate the literature references at the end of this paper. 
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2. Derivatives of the Gruneisen Function 

Defining x= 0/T, eq (2) becomes 

° W -x 4 Jo (e«-l)(l-«-t) 
Differentiating eq (3) and simplifying the result, one obtains 

r5 d£ 



j- % [G°(*)] = Gi(*)=-~ 



_4 r* rj\ 

x 4 )o (*<-!)( 



+ 



(l-e-«) J % L (l-e- r ) 2 



(3) 



(4) 



The Einstein function is defined as 



EHx) 



{\-e-*y 



(5) 



Substitution of eqs (3) and (5) in eq (4) yields 



G*(x)=-- [G°(x)-E°(x)]. 
x 



(6) 



For convenience, the variable x in the functions G°{x), E {) (x), and in their derivatives G 1 (x) , 
E 1 (x), etc. will be omitted and will be written simply as G°, E°, G 1 , £\ etc., respectively. One may 
note that the superscript numerals on both G and £ indicate the order of differentiation. 

The following is a list of the first five derivatives of the Gruneisen function obtained by the suc- 
cessive differentiation of eq (3) expressed in terms of the lower order derivatives of the Gruneisen 
and Einstein functions. 



4 4 

Gi = --G° + -E° 
x x 



G 2 =--G*+-E* 

X X 



6 4 

X X 



(7) 
(8) 
(9) 



7 4 

-G 3 + -£ 3 



G* = --G 4 + z £ 4 

* x 



(10) 

(ii) 



The recursion formula for the above is 



{ Z±l) G n-l + ± E *. 

X X 



(12) 



Where, n is the order of differentiation. It may be seen that the nth derivative of the Gruneisen 
function is given in terms of its (/? — l)th derivative. In order to be able to express the nth deriva- 
tive of the Gruneisen function in terms of the function itself, successive substitutions for the lower 
order derivatives are made. 



176 



The following is a list of the first five derivatives of the Griineisen function in terms of the 
Griineisen function and the derivatives of the Einstein function. 



G l = --(G°-E°) 



20 4 

G 2 = ^i (&-&)+- E* 



L20 24 4 



(6 



E l +-E 2 



X q X* x z X 






224 32 4 

^E 2 -^E* + -E* 



The recursion formula for the above is 



G n 



$>nC, n sqq_£q\ + SnzC ja ^ + S^C^ gg^ ^ | SniCni £j_, j + 



•5 nh'> /?/< 



(13) 

(14) 
(15) 
(16) 

(17) 



£"»- 



(18) 



where, C ,„■ and ,S„; indicate the absolute value and the sign of the coefficients, respectively, and are 
expressed by the following relations. 



C„ = 4 

C„, = [i>»] 

Ci=(n + 3)C»-i., L»^2] 

6'„i = (/? — /') C„_i,i + C B -i,i-i [n 25 / > 2] 

S«,= (-1)" 

S„, = (-l)" +i [n 55(^2] 



(19) 



The numerical values of the product SniCni required for the first seven derivatives of the 
Griineisen function are given in table 1. 

It can be seen from eq (18) that the relation for the derivative of the Griineisen function con- 
tains the derivatives of the Einstein function. The latter are obtained in the next section. 

3. Derivatives of the Einstein Function 

The Einstein function is defined by eq (5). For mathematical convenience, it may be written as 



£0 = 



1-e- 



y 



(20) 
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Table 1. List of the coefficients S ni C ni that appear in eq (18) 




-4 




















20 


4 

















-120 


-24 


4 














840 


168 


-28 


4 











-6,720 


-1,344 


224 


-32 


4 








60,480 


12,096 


-2,016 


288 


-36 


4 





604,800 


-120,960 


20,160 


-2,880 


360 


-40 


4 



where, 



y~- 



\- e -x 

Differentiation of eq (20) and simplification yields 



1 

' e x -\ 



^f- = E 1 = 2D E° 
ax 



where, 



The second derivative of E° becomes 



dE 



D° 



H- 



dx 



where, 



D 1 



dD° 
dx 



-\ + y+y\ 



(21) 



(22) 



(23) 



(24) 



(25) 



It should be noted that the superscript numeral of D indicates the order of differentiation 
similar to those used for G and E. It may be seen that the derivatives of the Einstein function con- 
tain the derivatives of the quantity y m for various values of m, where m indicates the power of y. 
The recursion formula for the derivatives of y" 1 is 



d_ 
dx 



( y m)=_ my m( y+ l)_ 



(26) 



The following is a list of the first five derivatives of the Einstein function obtained by the 
successive differentiation of eq (20). 



£ 1 = 2£>°£ 

E 2 =2D°E i + 2D i E° 

£ 3 = 2D°E 2 + WE 1 + 2D 2 E° 



(27) 
(28) 
(29) 
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E* = 2D°E 3 4- 6D <E Z + 6D 2 E' + 2D*E° 
E*=2D°E*+8D 1 E*+ 12D*E* t 8D*E 1 +2D*E° 

The recursion formula for the above is 

E"=K in D E'>~ l + K,r>D<E" -*+ K„J)*E»-*+ . . . + K lli D'"- l E"~ ! + . 

where, n indicates the order of differentiation. 

K„i are constants and can be obtained from the following relations 

K„\ = 2 

K ni = [i>n] 

Kni = K n -i,i + K n -i,i-i [n^i^2] 



+ K„ ll D"- i E° 



(30) 
(31) 

(32) 



(33) 



D m are coefficients which are functions of x. The superscript m indicates the order of differen- 
tiation. They are expressed by the following equations. 



»«=-\+\~y 



(34) 



D m = (— 1)^+1 
where, 



777 4- K mX y + Knny 2 + K m: \y :i + . . . + k,,,^' + . . . + K m , m+ iy m t 3 



[ m 5= 1 ] 



K m \ = K12 = 1 

* mj = [/> (i» + l)] 

Kmj= j(K m - Uj ) + 0'- l)(K m - W -i) [/ ^ 2, m ^ 2, J ^ (lH+1)]. 



(35) 



(36) 



The coefficients K n \ and k„,j can be represented as matrix elements. Their numerical values, 
required for the first seven derivatives of the Einstein function, are given in tables 2 and 3, 
respectively. 

Table 2. List of the coefficients K ni that appear in eq (32) 











i — * 








K„i 


1 


2 


3 


4 


5 


6 


7 


1 


2 




















2 


2 


2 

















n 3 


2 


4 


2 














1 4 


2 


6 


6 


2 











5 


2 


8 


12 


8 


2 








6 


2 


10 


20 


20 


10 


2 





7 


2 


12 


30 


40 


30 


12 


2 





Table 3. 


List of the coefficients K mi that 


appear in eq 


(35) 












J -* 










Knij 


1 2 


3 4 


5 


6 




7 


8 


1 


1 1 



















2 


1 3 


2 
















m 3 


1 7 


12 6 
















1 4 


1 15 


50 60 


24 













5 


1 31 


180 390 


360 


120 










6 


1 63 


602 2,100 


3,360 


2,520 




720 





7 


1 127 


1,932 10,206 


25,200 


31,920 


20,160 


5,040 
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388-513 O - 70 -4 



Calculating the derivatives of the Einstein function using eq (32) and substituting these in 
either eq (12) or eq (18) one can obtain the derivatives of the Griineisen function. Since the nth 
derivative of the Griineisen function contains all lower derivatives (up to n — 1) of the Einstein 
function, it is sufficient to derive a recursion formula for the derivatives of the Einstein function 
in terms of its lower derivatives. However, in the case of the Griineisen function, a recursion 
formula for its derivatives was obtained both in terms of lower derivatives of G and in terms of G 
itself. 

4. Computation of the Derivatives of the Einstein and Griineisen Functions 

Based on the formulas developed in the previous sections, the first five derivatives of the 
Einstein and Griineisen functions are computed for the range 1 < x < 30 using a computer. The 
selection of the limits is that of a practical one. Above the Debye temperature (x = 1) both functions 
are relatively insensitive to temperature variations and approach unity as temperature increases. 
Below T= 0/S0(x = 30) the values of the functions become small and also the integral in the 
Griineisen function attains a constant value. 

The Einstein function and its first five derivatives are computed using eqs (20), (27), (28), (29), 
(30), and (31), respectively in conjunction with eqs (34), (35), and (36). The results are given in 
table 4. The values obtained for the Einstein function are in exact agreement with those reported by 
Hilsenrath and Ziegler [3]. 









Table 4. The Einstein function and its first five derivatives* 








y 


E° 




E l E 2 E s E* 




£ 5 




A 


Coef 


Exp 


Coef Exp Coef Exp Coef Exp Coef 


Exp 


Coef 


Exp 



1 


9.20674 


-1 


- 1.50948 


— 1 


- 1.21319 


- 1 


8.19752 


-2 


4.93166 


-2 


- 8.53277 


-2 


2 


7.24062 


-1 


- 2.26657 


— 1 


- 2.89466 


-2 


8.94460 


-2 


- 2.89868 


-2 


- 5.50503 


-2 


3 


4.96269 


-1 


-2.17428 


— 1 


3.97080 


-2 


4.43374 


-2 


- 5.06061 


-2 


7.34212 


-3 


4 


3.04087 


-1 


- 1.63391 


— l 


6.13399 


-2 


2.48308 


-3 


- 3.00279 


-2 


2.62514 


-2 


5 


1.70742 


- 1 


- 1.04762 


— 1 


5.29512 


-2 


- 1.54893 


-2 


- 7.59110 


-3 


1.66883 


-2 


6 


8.96791 


-2 


- 6.02318 


-2 


3.59185 


-2 


- 1.68202 


-2 


3.02914 


-3 


5.38369 


-3 


7 


4.47638 


-2 


- 3.20559 


-2 


2.12103 


-2 


- 1.22491 


-2 


5.20388 


-3 


- 1.14766 


-4 


8 


2.14840 


-2 


- 1.61274 


-2 


1.14494 


-2 


- 7.45505 


-3 


4.15200 


-3 


- 1.55181 


-3 


9 


9.99866 


-3 


- 7.77921 


-3 


5.80800 


-3 


- 4.08606 


-3 


2.61503 


-3 


-1.39754 


-3 


10 


4.54041 


-3 


- 3.63274 


-3 


2.81612 


-3 


- 2.09076 


-3 


1.45695 


-3 


-9.15259 


-4 


11 


2.02097 


-3 


- 1.65359 


-3 


1.31966 


-3 


- 1.01920 


-3 


7.52291 


-4 


-5.19021 


-4 


12 


8.84777 


-4 


- 7.37325 


-4 


6.02169 


-4 


-4.79315 


-4 


3.68772 


-4 


-2.70560 


-4 


13 


3.81997 


-4 


- 3.23230 


-4 


2.68985 


-4 


-2.19263 


-4 


1.74065 


-4 


- 1.33395 


-4 


14 


1.62980 


-4 


- 1.39697 


-4 


1.18078 


-4 


-9.81221 


-5 


7.98299 


-5 


- 6.32020 


-5 


15 


6.88281 


-5 


- 5.96510 


-5 


5.10858 


-5 


- 4.31325 


-5 


3.57911 


-5 


-2.90616 


-5 


16 


2.88090 


-5 


- 2.52079 


-5 


2.18318 


-5 


- 1.86809 


-5 


1.57550 


-5 


- 1.30542 


-5 


17 


1.19644 


-5 


- 1.05568 


-5 


9.23206 


-6 


- 7.99009 


-6 


6.83091 


-6 


-5.75454 


-6 


18 


4.93451 


-6 


- 4.38623 


-6 


3.86842 


-6 


- 3.38106 


-6 


2.92416 


-6 


- 2.49772 


-6 


19 


2.02261 


-6 


- 1.80970 


-6 


1.60800 


-6 


- 1.41751 


-6 


1.23822 


-6 


- 1.07013 


-6 


20 


8.24461 


-7 


- 7.42015 


-7 


6.63691 


-7 


- 5.89490 


-7 


5.19411 


-7 


-4.53454 


-7 


21 


3.34391 


-7 


- 3.02544 


-7 


2.72214 


-7 


- 2.43400 


-7 


2.16103 


-7 


- 1.90322 


-7 


22 


1.35010 


-7 


- 1.22737 


-7 


1.11021 


-7 


- 9.98630 


-8 


8.92630 


-8 


- 7.92209 


-8 


23 


5.42853 


-8 


- 4.95649 


-8 


4.50497 


-8 


- 4.07397 


-8 


3.66349 


-8 


- 3.27354 


-8 


24 


2.17448 


-8 


- 1.99327 


-8 


1.81961 


-8 


- 1.65351 


-8 


1.49495 


-8 


-1.34395 


-8 


25 


8.67997 


-9 


- 7.98557 


-9 


7.31895 


-9 


- 6.68010 


-9 


6.06903 


-9 


- 5.48574 


-9 


26 


3.45374 


-9 


-3.18807 


-9 


2.93262 


-9 


- 2.68738 


-9 


2.45236 


-9 


- 2.22756 


-9 


27 


1.37018 


-9 


- 1.26868 


-9 


1.17095 


-9 


- 1.07697 


-9 


9.86753 


- 10 


-9.00294 


- 10 


28 


5.42089 


-10 


- 5.03368 


-10 


4.66031 


- 10 


- 4.30076 


- 10 


3.95504 


- 10 


-3.62315 


-10 


29 


2.13922 


-10 


- 1.99169 


-10 


1.84925 


-10 


-1.71189 


- 10 


1.57962 


- 10 


- 1.45243 


- 10 


30 


8.42186 


- 11 


- 7.86040 


-11 


7.31766 


- 11 


- 6.79363 


- 11 


6.28832 


- 11 


-5.80173 


- 11 



* The superscript numerals indicate the order of differentiation. 



The Griineisen function and its first five derivatives are computed using eqs (3), (13), (14), 
(15), (16), and (17), respectively. The results are given in table 5. The Griineisen function expressed 
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TABLE 5. The Griineisen function and its first five derivatives* 














jv 




G° 




G l G 2 G :i 




G 4 




G* 




A 


Coef 


Exp 


Coef 


Exp 


Coef Exp Coef Exp Coef 


Exp 


Coef 


Exp 


Coef 


Exp 



GO 


1 


2.36616 


-1 


9.46463 


— 1 


-1.03156 


— 1 


-8.80120 


-2 


4.27961 


-2 


2.83281 


-2 


-2.93585 


-2 


O 


2 


3.22929 





8.07321 


— 1 


-1.66518 


— 1 


-3.70190 


-2 


5.31638 


-2 


-7.18128 


-3 


-2.92485 


-2 


1 


3 


1.27705 


+ 1 


6.30643 


— l 


-1.79165 


— 1 


8.70489 


-3 


3.55342 


-2 


-2.37967 


-2 


-4.01706 


-3 


o 


4 


2.94882 


+ 1 


4.60754 


— 1 


-1.56667 


— 1 


3.24428 


-2 


1.26758 


-2 


-1.96995 


-2 


9.37115 


-3 


en 


5 


5.02629 


+ 1 


3.21682 


— 1 


-1.20752 


— 1 


3.69424 


-2 


-1.96992 


-3 


-9.63355 


-3 


9.34080 


-3 




6 


7.08729 


+ 1 


2.18744 


— i 


-8.60433 


-2 


3.15482 


-2 


-7.60252 


-3 


-2.34386 


-3 


5.14457 


-3 




7 


8.83248 


+ 1 


1.47147 


— 1 


-5.85047 


-2 


2.34714 


-2 


-7.99817 


-3 


9.98689 


-4 


1.83229 


-3 




8 


1.01477 


+ 2 


9.90984 


-2 


-3.88072 


-2 


1.61908 


-2 


-6.41840 


-3 


1.88858 


-3 


1.87425 


-4 




9 


1.10546 


+ 2 


6.73958 


-2 


-2.55098 


-2 


1.07147 


-2 


-4.56180 


-3 


1.73204 


-3 


-3.77359 


-4 




10 


1.16379 


+ 2 


4.65517 


-2 


-1.68045 


-2 


6.94916 


-3 


-3.04305 


-3 


1.29383 


-3 


-4.52284 


-4 




11 


1.19927 


+ 2 


3.27647 


-2 


-1.11795 


-2 


4.48030 


-3 


-1.96393 


-3 


8.79152 


-4 


-3.65823 


-4 




12 


1.21988 


+ 2 


2.35317 


-2 


-7.54897 


-3 


2.89963 


-3 


-1.24909 


-3 


5.68866 


-4 


-2.56320 


-4 




13 


1.23141 


+ 2 


1.72460 


-2 


-5.18892 


-3 


1.89628 


-3 


-7.92444 


-4 


3.59235 


-4 


-1.67509 


-4 




•— ' 14 

22 Is 


1.23765 


+ 2 


1.28868 


-2 


-3.63538 


-3 


1.25844 


-3 


-5.05593 


-4 


2.24762 


-4 


-1.05627 


-4 




1.24093 


+ 2 


9.804S9 


-3 


-2.59628 


-3 


8.49521 


-4 


-3.26185 


-4 


1.40718 


-4 


-6.55052 


-5 




16 


1.24262 


+ 2 


7.58434 


-3 


-1.88888 


-3 


5.83974 


-4 


-2.13532 


-4 


8.87501 


-5 


-4.04363 


-5 




17 


1.24347 


+ 2 


5.95523 


-3 


-1.39842 


-3 


4.08815 


-4 


-1.42115 


-4 


5.66380 


-5 


-2.50459 


-5 




18 


1.24388 


+ 2 


4.73969 


-3 


-1.05217 


-3 


2.91294 


-4 


-9.62384 


-5 


3.66747 


-5 


-1.56501 


-5 




19 


1.24409 


+ 2 


3.81853 


-3 


-8.03475 


-4 


2.11060 


-4 


-6.63120 


-5 


2.41323 


-5 


-9.90029 


-6 




20 


1.24418 


+ 2 


3.11045 


-3 


-6.21925 


-4 


1.55333 


-4 


-4.64671 


-5 


1.61456 


-5 


-6.35436 


-6 




21 


1.24422 


+ 2 


2.55907 


-3 


-4.87378 


-4 


1.15985 


-4 


-3.30867 


-5 


1.09825 


-5 


-4.14266 


-6 




22 


1.24424 


+ 2 


2.12459 


-3 


-3.86265 


-4 


8.77651 


-5 


-2.39157 


-5 


7.59140 


-6 


-2.74428 


-6 




23 


1.24425 


+ 2 


1.77851 


-3 


-3.09297 


-4 


6.72298 


-5 


-1.75304 


-5 


5.32824 


-6 


-1.84693 


-6 




24 


1.24425 


+ 2 


1.50011 


-3 


-2.50015 


-4 


5.20831 


-5 


-1.30177 


-5 


3.79408 


-6 


-1.26220 


-6 




25 


1.24425 


-f 2 


1.27411 


-3 


-2.03856 


-4 


4.07700 


-5 


-9.78362 


-6 


2.73834 


-6 


-8.75299 


-7 




26 


1.24425 


+ 2 


1.08912 


-3 


-1.67556 


-4 


3.22219 


-5 


-7.43537 


-6 


2.00142 


-6 


-6.15443 


-7 




27 


1.24425 


+ 2 


9.36511 


-4 


-1.38742 


-4 


2.56928 


-5 


-5.70934 


-6 


1.48004 


-6 


-4.38384 


-7 




28 


1.24425 


+ 2 


8.09722 


-4 


-1.15674 


-4 


2.06561 


-5 


-4.42624 


-6 


1.10650 


-6 


-3.16086 


-7 




29 


1.24425 


+ 2 


7.03681 


-4 


-9.70594 


-5 


1.67344 


-5 


-3.46225 


-6 


8.35693 


-7 


-2.30514 


-7 




30 


1.24425 


+ 2 


6.14445 


-4 


-8.19260 


-5 


1.36543 


-5 


-2.73085 


-6 


6.37190 


-7 


-1.69909 


-7 



*The superscript numerals indicate the order of differentiation. 
* Forx = oc, J,= 124.4313 [2]. 



by eq (3) may be written as 



^here 



G W=(^)A (37) 

J " I (e«-l)(l -«-<)" (38) 

There is no known closed form for the above integration, therefore numerical integration techniques 
are employed. Various values for the increment in numerical integration are tried to determine 
the convergence of the results. Values 0.002 and 0.0002 yielded a difference of 0.00042 percent 
for x = 1, and 0.0048 percent for x = 30 in y 5 . In the literature, Rogers and Powell [4] have computed 
the integral J 5 utilizing series expansion approximation. Their results correspond exactly to those 
obtained using a value of 0.002 for the increment. Since the accuracy of the approximation of nu- 
merical integration increases with decreasing increment size, a value of 0.0002 is used in the 
final computations. Thus, it is expected that the present results of J$ and hence the Griineisen 
function are more accurate than those given by Rogers and Powell. 

It may be noted that the only numerical method employed in the present computations is 
that for the evaluation of 7 5 , the integral appearing in the Griineisen function. All other computa- 
tions including those for the first five derivatives of the Griineisen and Einstein functions represent 
mathematically closed forms and thus are not subject to errors common to numerical differentiation. 

5. Summary 

The Griineisen and Einstein functions appear in the mathematical expressions describing a 
number of thermodynamic and transport properties of matter. In order to improve the interpreta- 
tion of experimental results and extend the application of the theories regarding these properties, 
a knowledge of the temperature derivatives of the above functions are needed. In the present 
study, mathematical expressions in closed form are derived for the derivatives of the Griineisen 
and Einstein functions. Also, the first five derivatives of both functions are computed and presented 
in tabular form. Since all the derivatives are obtained from mathematical relations of closed 
form, no errors inherent to numerical differentiation are introduced. 
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